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This  investigation  is  an  attempt  to  verify  the  results  of  a  theoretical  model, 
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Poisson’s  Ratio  was  calculated  from  the  shear  and  Young’s  moduli  measured  from  a 
cylindrical  sample  composed  of  300  |im  glass  beads.  This  was  used  to  calculate  the  bulk 
moduli  for  the  porous  spheres;  the  shear  moduli  h^d  been  previously  measured  by  LT. 
Huskey.  These  and  other  material  properties  were  used  as  input  to  the  theoretical  model 
developed  by  Kargl  and  Lim. 

The  experimental  data  were  compared  to  the  theoretical  values.  Reasonably  good 
agreement  between  the  measured  and  predicted  scattering  was  obtained  for  the  aluminum 
and  100  |xm  spheres.  The  measured  scattering  from  the  500  |im  sphere  was  in  poor 
agreement  with  the  theoretical  predictions. 
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I.  INTRODUCTION 


A.  BACKGROUND 

Because  of  the  increased  emphasis  on  littoral  warfare  and  the  threat  of  the  use  of 
buried  mines  in  these  areas,  there  is  increased  interest  in  the  scattering  of  sound  in  fluid- 
saturated  porous  media. 

Maurice  Biot  (Biot,  1956a,  b)  developed  a  general  theory  for  the  propagation  of 
elastic  waves  in  a  fluid-saturated  porous  media.  Biot’s  theory  has  been  applied  to  the 
scattering  of  elastic  waves  from  a  saturated  porous  sphere  in  a  saturated  porous  host 
(Kargl  and  Lim,  1993). 

In  1993,  LT.  Huskey  performed  experiments  in  an  attempt  to  verify  Kargl  and 
Lim  s  model.  The  results  of  his  experiment  were  good  when  10%  frame  damping  was 
included  in  the  calculations.  Other  calculations  were  less  satisfactory.  (Huskey,  1993) 

B.  OBJECTIVES 

Kargl  and  Lim  developed  a  numerical  model  to  compute  the  scattering  of  sound 
from  a  saturated  porous  sphere.  This  research  is  an  attempt  to  experimentally  measure  the 
scattering  of  sound  from  a  saturated  porous  sphere  and  compare  the  results  with  the 
numerical  values  obtained  from  Kargl  and  Lim’s  model  (Kargl  and  Lim,  1993). 

C.  EXPERIMENT  OVERVIEW 

As  stated  above,  the  purpose  of  this  research  is  to  measure  the  scattering  of  sound 
from  a  fluid  saturated  porous  sphere.  The  spheres  which  were  used  in  this  research  were 
composed  of  borosilicate  glass  beads.  The  beads  had  a  mean  diameter  of  100  |im  and 
500  pm.  The  beads  were  coated  with  a  heat  curing  epoxy  powder  and  then  poured  into 
cylindrical  molds.  These  molds  were  then  heated  to  cure  the  epoxy.  The  resulting 
cylinders  were  ground  into  spheres  with  a  diameter  of  about  6.8  centimeters.  Cylindrical 
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rods  composed  of  100  pm  and  500  pm  glass  beads  were  made  at  the  same  time  as  the 
spheres.  These  rods  were  used  by  LT.  Huskey  to  measure  the  permeability,  porosity,  and 
shear  modulus  of  the  porous  matrices.  These  values  are  summarized  in  Table  1  (Huskey, 
1993). 


Sample 

100  pm 

500  pm 

Permeability,  k 

6.53x1 0-'^m^ 

5.74x10-”  m® 

Porosity,  P 

0.309 

0.321 

Shear  Modulus,  G 

2.81x10®  Pa 

2.72x10®  Pa 

Table  1.  Material  measurements  taken  by  LT.  Huskey. 


Scattering  measurements  taken  for  these  spheres  were  compared  to  theoretical 
values  determined  by  the  model  developed  by  Kargl  and  Lim.  This  model  requires  eleven 
material  parameters  to  determine  the  scattered  amplitude.  The  scattering  solid  provides 
five  of  these  material  parameters:  mass  density,  bulk  and  shear  moduli  of  the  solid,  and 
bulk  and  shear  moduli  of  the  porous  lattice.  Three  material  properties  are  determined  by 
the  structure  of  the  porous  lattice;  tortuosity,  permeability,  and  structural  factor.  The  final 
three  properties  are  of  the  bulk  fluid;  mass  density,  and  the  bulk  and  shear  moduli.  (Kargl 
and  Lim).  The  program  provided  by  Kargl  used  the  bulk  fluid’s  bulk  modulus  and 
viscosity  to  determine  complex  values  for  the  bulk  and  shear  moduli. 

In  this  experiment  all  the  material  properties  were  accurately  known  with  the 
exception  of  the  bulk  moduli  of  the  spheres.  This  could  not  be  measured  by  LT.  Huskey 
at  the  time  of  his  research  (Huskey,  1993)  due  to  the  small  aspect  ratio  (length  to 
diameter)  of  the  cylindrical  samples.  For  this  experiment  a  new  cylinder  of  300  pm 
diameter  porous  glass  beads  with  greater  aspect  ratio  was  obtained  so  that  the  bulk  and 
shear  moduli  could  be  determined.  From  these  moduli  the  Poisson’s  Ratio  can  be 
determined  which  can  then  be  used  along  with  the  shear  moduli  of  the  100  and  500  pm 
samples  to  determine  their  bulk  moduli.  Once  these  parameters  are  determined  they  can 
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be  used  as  input,  along  with  the  other  known  material  properties,  to  Kargl  and  Lim’s 
theoretical  model  and  the  results  compared  with  the  experimentally  measured  values. 
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n.  DETERMINATION  OF  ELASTIC  MODULI 


A.  RESONANT  ACOUSTIC  METHOD 

The  elastic  moduli  of  a  cylindrical  sample  can  be  obtained  using  a  resonant 
acoustic  method  (Garrett,  1990).  This  method  employs  a  transducer  bonded  to  each  end 
of  a  cylindrical  sample  which  should  have  a  length  to  diameter  ratio  »1.  The  transducers 
are  made  from  coiled  magnet  wire  and  are  attached  with  epoxy.  These  are  used  to  set  up 
flexural,  torsional,  and  longitudinal  standing  waves  in  the  samples.  The  bars  are 
positioned  so  that  the  transducers  attached  to  the  ends  are  centered  between  the  pole  faces 
of  strong  magnets.  Depending  on  the  orientation  of  the  transducers  to  the  pole  pieces,  any 
of  the  three  vibrational  modes  can  be  selectively  excited.  Figure  1,  from  (Garret,  1990), 
illustrates  the  position  of  the  magnet  faces  for  the  torsional  and  longitudinal  modes.  The 
fundamental  frequency  of  each  mode  can  be  determined  based  on  the  length  and 
boundary  conditions  of  the  sample.  These  frequencies  can  then  be  used  to  determine  the 
moduli  of  the  sample. 

In  the  experiments  conducted  by  LT.  Huskey  only  the  shear  modulus  could  be 
measured  due  to  the  length  to  diameter  ratio  of  the  samples  (about  3:1)  which  intensified 
electrical  cross-talk  between  the  attached  transducers. 


Figure  1.  Transducer  placement  for  the  torsional  and  longitudinal  modes. 
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B.  ELASTIC  MODULI  MEASUREMENT 


A  sample  made  of  300  |u,m  glass  beads  with  a  length  of  19.7  centimeters  and 
diameter  of  2.6  centimeters  (giving  a  length  to  diameter  ratio  of  7.6:1)  was  obtained  and 
the  resonant  acoustic  method  was  used  to  determine  both  the  shear  and  bulk  moduli  of 
the  sample.  Approximately  2  meters  of  #32  magnet  wire  were  attached  to  each  end  of  the 
sample  with  epoxy  resulting  in  15  turns  of  wire  in  each  transducer.  Care  was  taken  to 
minimize  the  amount  of  weight  added  by  the  epoxy.  The  sample  weighed  153.52  grams 
and  the  transducers  added  4.38  grams.  The  mass  loading  of  the  transducers  was 
significant  enough  that  an  effective  length  of  the  bar  had  to  be  calculated  using  the 
following  equations  (Garret,  1990): 


=  hi  1  +  (longimdinal  mode) 


MJ 

=  £|  1  +  I  (torsional  mode) 


where  L  is  the  actual  length  of  the  sample,  m  is  the  transducer  mass,  and  M  is  the  mass  of 
the  sample. 

The  sample  was  positioned  so  that  each  transducer  was  centered  between  the  pole 
pieces  of  the  magnets.  One  set  of  magnets  were  rotated  90°  to  minimize  electrical  cross¬ 
talk  between  the  two  transducers  (one  transducer  was  also  offset  90°  relative  to  the  other 
when  it  was  attached).  A  Hewlett-Packard  35665A  Dynamic  Signal  Analyzer  was  used  to 
drive  one  transducer  through  a  Hewlett-Packard  467A  Power  Amplifier.  The  signal  from 
the  other  transducer  was  then  input  to  the  Dynamic  Signal  Analyzer.  Both  signals  were 
monitored  on  an  oscilloscope  to  ensure  that  no  distortion  was  occurring.  The  output  from 
the  Dynamic  Signal  Analyzer  was  then  swept  upward  in  frequency  from  100  Hz  to  20 
kHz.  The  resonances  were  noted  and  the  frequency  bracket  around  each  was  tightened  to 
allow  for  more  accurate  measurement  of  the  resonance. 
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_ Modal  Frequency  Summary  _ 

Mode  Number,  Frequency  (Hz),  Normalized  Frequency  (Hz), 

n _ fn _ tn/n _ 

Torsional _ _ _ 

1  _ 2505 _ 2505 _ 

2  _ 5175 _ 2588 _ 

3  7830  2610 


Table  2.  Resonant  frequency  measurements. 


The  first  four  modes  of  both  the  longitudinal  and  torsional  modes  were  measured 
and  averaged  together.  The  results  are  shown  in  Table  2.  The  Young’s  modulus  (E)  of  the 
sample  was  determined  using  the  following  equation  (from  Garrett): 

where  p  is  the  mass  density,  n  is  the  mode  number,  and//  is  the  frequency  of  the  nth 
longitudinal  mode.  The  shear  modulus  (G)  was  similarly  calculated  using  (from  Garrett) 

G  =  /nf 

where//  is  the  frequency  of  the  nth  torsional  mode. 

The  calculated  Young’s  modulus  was  3.83x10^  Pa  and  the  shear  modulus  was 
1.68x10^  Pa.  These  values  are  approximately  1/20  of  the  borosilicate  glass  used  to  make 

the  glass  beads.  These  were  used  to  calculate  the  Poisson’s  Ratio  (v)  by  combining 

^  3K-2G 
2{3K  +  G) 


and 

K-  ^ 

3(1 -2v) 

to  give 

< 

II 

1 
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Shear  Modulus,  G, 

(from  Huskey) 

100  urn 

500  pm 

2.81x10®  Pa 

2.72x10®  Pa 

"Bulk  Modulus,  K, 
assuming  v  =  0.14. 

2.96x10®  Pa 

2.87x10®  Pa 

Table  3.  Material  properties  of  the  samples. 


where  K  is  the  bulk  modulus.  The  resulting  Poisson’s  Ratio  was  0. 14.  This  is  an 
unusually  small  Poisson’s  Ratio  and  is  possibly  due  to  the  epoxy  used  to  bond  the  glass 
beads  together.  The  epoxy  may  be  stretching  between  the  beads  allowing  the  cylinder  to 

elongate  without  much  lateral  constriction. 

Since  the  300  [im  rod  was  made  at  a  different  time  than  the  100  |im  and  500  pm 
spheres,  the  shear  moduli  measured  by  Huskey  and  the  calculated  Poisson’s  Ratio  were 
used  to  calculate  new  bulk  moduli  for  both  spheres.  These  values  are  listed  in  Table  3  and 
were  used  as  inputs  to  Kargl’s  program. 


8 


in.  BACKSCATTER  MEASUREMENTS 


A.  MEASUREMENT  OBJECTIVE 

The  objective  of  the  backscatter  measurement  was  to  separate  the  scattered 
acoustic  pressure  from  the  multipath  interference.  Multipath  interference  was  caused  by 
surface  reflections  from  the  transmitter’s  side  lobes  and  from  piping  used  for  filtration  of 
the  tank’s  water. 

B.  EXPERIMENT  SETUP 

The  backscatter  measurements  were  taken  in  a  water-filled  tank,  in  Spanagel  Hall 
Room  025,  measuring  7.3  meters  in  length,  1.6  meters  in  width  and  2.0  meters  in  depth. 
The  walls  and  bottom  of  the  tank  are  covered  with  anechoic  tile.  Figure  2  shows  the  tank 
setup  for  the  backscatter  measurements.  All  components  were  aligned  along  the 
centerline  of  the  tank  as  viewed  from  the  top  of  the  tank. 

A  type  F33  general-purpose  directional  transducer  was  used  as  the  projector.  This 
transducer  is  shown  in  Figure  3  from  the  USED  Transducer  Catalog,  April  1991.  Its 


2.0  m 


■2.38  m 


1.1  m 


Type  F-33 
Transducer 


-0.75  m-^ 


B&K8103  Sphere 
Hydrophone 


'7.3  m- 


Figure  2.  Tank  setup  for  scattering  measurements. 
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design  consists  of  a  small  inner  array  made  of  12  PZT  disks  with  a  frequency  range  of 
15-150  kHz.  The  larger  outer  array  consists  of  64  PZT  squares  with  a  frequency  range  of 
1-50  kHz  (USRD  Transducer  Catalog).  The  arrays  can  be  used  individually  or  wired  in 
parallel.  In  this  experiment  the  arrays  were  connected  in  parallel  for  added  directionality. 

To  ensure  accurate  measurements,  the  scatterer  must  be  placed  in  the  far  field  of 
the  transmitter.  The  following  equation  was  used  to  determine  the  limiting  distance  to  the 
far  field  (r^J: 

min  ji 

4  A 

where  d  is  the  dimension  of  the  transmitter  and  A  is  the  wavelength  of  sound  (Kinsler  et 
al.).  The  lowest  frequency  used  in  this  experiment  was  30  kHz,  giving  a  maximum 
distance  to  the  far  field  of  less  than  0.25  meters. 
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The  receiver  used  was  a  Bruel  and  Kjaer  type  8103  hydrophone.  Figure  4  shows 
an  illustration  of  the  receiver  from  its  calibration  chart.  It  has  a  frequency  range  of  0. 1  to 
180  kHz.  The  receiver  was  suspended  from  the  arm  of  a  protractor-like  device  illustrated 
in  Figure  5.  The  arm  can  be  rotated  in  five  degree  increments  and  locked  into  place  with  a 
lock  pin.  The  180°  position  corresponds  to  sound  being  scattered  from  the  sphere  directly 
back  at  the  transmitter.  The  receiver  could  be  moved  radially,  in  increments  of  5 
centimeters,  between  40  to  75  centimeters  from  the  center  of  the  protractor.  A  small 
clamp  was  attached  to  the  cable  of  the  hydrophone  allowing  the  depth  to  be  adjusted.  The 
clamp  rested  in  any  one  of  several  beveled  holes  in  the  protractor’s  arm.  A  weight  was 


Units  are  mm 

Figure  4.  Bruel  &  Kjaer  Type  8103  hydrophone. 


suspended  from  the  bottom  of  the  receiver  to  ensure  that  it  hung  straight  down.  The 
weight  was  suspended  a  few  centimeters  from  the  bottom  of  the  tank. 

The  target  spheres  were  suspended  in  a  fine  net  by  a  string.  The  string  passed 
through  a  hole  drilled  in  the  bolt  which  was  the  pivot  point  of  the  protractor’s  arm.  This 
allowed  the  depth  of  the  sphere  to  be  adjusted  to  correspond  to  the  depth  of  the  center  of 
the  transmitter.  The  spheres  were  degassed  to  ensure  that  no  air  bubbles  were  trapped  in 
the  spheres  during  the  experimental  measurements.  This  was  done  by  placing  them  in  a 


beaker  of  water  and  then  placing  the  beaker  under  a  bell  jar.  A  vacuum  pump  was  then 
used  to  evacuate  the  bell  jar  until  the  water  began  to  boil.  At  this  point  the  hose  to  the  bell 
jar  was  clamped  and  the  vacuum  pump  was  turned  off.  Periodically  the  beaker  and  sphere 
were  agitated  to  dislodge  any  bubbles  adhering  to  the  surface.  The  vacuum  was  held  for 
approximately  24  hours  to  ensure  no  bubbles  remained  trapped  in  the  spheres.  Next  the 
hose  clamp  to  the  bell  jar  was  removed  allowing  the  pressure  to  return  to  normal.  The 
bell  jar  was  carefully  removed  and  the  sphere  transported  to  the  water  tank  in  the  beaker 
full  of  water.  The  beaker  was  submerged  in  the  water  tank  and  the  sphere  removed  thus 
keeping  it  submerged  at  all  times. 


Figure  5.  Protractor  device  used  to  position  receiver. 
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Figure  6.  Electronic  equipment  setup. 


Figure  6  is  a  block  diagram  of  the  electronics  used  for  this  experiment.  A  1  Volt 
peak  to  peak  sine  wave  was  generated  with  the  Hewlett-Packard  3314A  Signal  Generator 
which  was  then  amplified  to  a  10  Volt  peak  to  peak  signal  by  a  Hewlett-Packard  467A 
Power  Amplifier.  This  signal  was  then  applied  to  the  Type  F33  transducer.  The  signal 
could  also  be  monitored  by  either  the  Nicolet  Pro  30  Digital  Oscilloscope  or  the 
DSS5020  Oscilloscope.  The  signal  received  by  the  Bruel  &  Kjaer  8103  hydrophone  was 
input  to  an  Ithaco  1201  Preamplifier.  The  output  from  the  preamplifier  was  then  analyzed 
by  the  Nicolet  Pro  30  Oscilloscope. 

C.  EXPERIMENT4L  PROCEDURE 

The  frequency  range  used  for  this  experiment  was  30  kHz  to  150  kHz, 
corresponding  to  a  range  of  ka  of  approximately  4  to  22,  referred  to  the  wave  number  in 
water. 
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1.  Monostatic  Measurements 

For  the  monostatic  measurements  the  sphere  was  hung  below  the  protractor.  The 
B&K  hydrophone  was  positioned  75  centimeters  from  the  sphere  along  the  arm  of  the 
protractor  which  was  positioned  at  180°.  The  signal  generator  was  setup  to  send  out  30 
kHz  bursts  at  a  rate  of  five  bursts  per  second.  Each  burst  consisted  of  a  30-cycle  sine 
wave. 

The  Nicolet  Pro  30  was  used  to  monitor  the  received  signal.  A  trigger  delay  was 
used  so  that  the  trace  on  the  oscilloscope  would  begin  just  before  the  arrival  of  the 
scattered  signal.  The  delay  was  determined  by  the  speed  of  sound  and  the  distance 
between  the  transmitter,  sphere  and  receiver.  The  sample  rate  was  set  so  that  the  pulse 
would  fill  as  much  of  the  oscilloscope’s  screen  as  possible  without  being  cutoff.  Careful 
attention  was  paid  to  ensure  that  the  sample  rate  remained  well  above  the  Nyquist 
frequency  of  the  received  signal.  Fifty  bursts  of  the  received  signal  were  averaged 
together  and  then  saved  to  floppy  disk.  Measurements  were  taken  between  30  and  150 
kHz  in  2  kHz  steps. 

At  this  point  the  sphere  was  carefully  removed  from  the  tank  to  a  bucket  of  water 
and  the  above  procedure  was  repeated  to  obtain  the  multipath  interference  background  in 
the  tank.  The  next  measurement  to  be  taken  was  the  incident  signal  on  the  sphere.  Due  to 
the  design  of  the  protractor,  the  B&K  hydrophone  could  not  be  positioned  in  the  same 
place  as  the  sphere  but  had  to  be  positioned  15  centimeters  behind  it.  This  displacement 
was  taken  into  account  in  the  calculations  by  using  a  1/r  signal  fall-off  in  the  far  field  of 
the  transmitter  (Kinsler  et  al.).  The  incident  signal  was  measured  at  each  frequency. 

Once  the  received  signal  and  multipath  interference  background  measurements 
had  been  taken,  the  saved  waveforms  were  subtracted  from  each  other,  leaving  only  the 
scattered  pulse.  The  subtraction  process  was  performed  using  the  Nicolet  Pro  30  Digital 
Oscilloscope.  Figure  7  shows  the  measured  signal,  multipath  interference  background, 
and  resulting  scattered  signal  from  the  subtraction  process. 
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Multipath  interference  background  signal. 


Figure  7.  Example  of  subtracting  the  measured  signal  from  the  multipath  interference 
background.  At  top  is  the  received  signal,  in  the  middle  is  the  multipath  interference,  and  bottom 
is  the  difference  between  the  measured  signal  and  multipath  interference  background. 
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Next  a  Fast  Fourier  Transform  (FFT)  was  taken  of  the  scattered  and  incident 
waveforms,  using  the  Nicolet  Pro  30,  to  determine  the  amplitude  of  the  desired 
transmitted  frequency  component.  The  start  and  stop  points  of  the  FFT  were  chosen  so 
that  they  were  inside  the  start  and  end  points  of  the  transmitted  pulse.  This  was  done  to 
avoid  distortion  caused  by  the  HP467A  Power  Amplifier  turning  on  and  off,  and  to  avoid 
any  ring  up  and  down  of  the  Type  F33  transmitter.  Care  was  taken  to  ensure  that  a  whole 
number  of  wavelengths  were  taken  and  that  the  start  and  stop  points  were  as  close  as 
possible  to  a  zero  crossing  to  avoid  leakage  into  adjacent  frequency  bins  when  the  FFT 
was  taken  (Hewlett-Packard,  Application  Note  243,  pp.  25-26).  No  windowing  function 
was  used  with  the  above  method.  The  results  of  these  measurements  are  listed  in 
Appendix  A.  The  method  of  calculation  of  the  tabulated  results  are  discussed  in  Chapter 
IV  Section  A,  Normalization  of  Experimental  Data. 

2.  Bistatic  Measurements 

Bistatic  measurements  were  conducted  in  a  manner  similar  to  the  monostatic 
measurements.  The  difference  was  that  instead  of  adjusting  the  frequency  between 
measurements,  the  angular  position  of  the  receiver  was  adjusted.  The  receiver  was  started 
in  the  90°  position  and  then  moved  in  5°  increments  until  the  270°  position  was  reached. 
The  sphere  was  then  removed  from  the  tank  and  the  background  measurements  were 
taken.  Next  the  receiver  was  moved  to  measure  the  incident  signal  level  on  the  sphere. 
After  these  measurements  were  made  the  resulting  waveforms  were  subtracted  and  FFTs 
taken  to  obtain  the  scattered  and  incident  signal  levels.  These  measurements  were 
performed  at  the  following  frequencies;  30,  60,  90,  120,  and  150  kHz.  Results  of  the 
measurements  are  listed  in  Appendix  B. 

For  the  aluminum  sphere  measurements  were  taken  only  at  30  kHz  and  150  kHz 
since  these  were  to  be  used  as  a  gauge  of  the  effectiveness  of  the  experimental  method. 
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IV.  DATA  ANALYSIS  AND  RESULTS 


A.  NORMALIZATION  OF  EXPERIMENTAL  DATA 


All  of  the  experimental  data  was  normalized  to  a  distance  of  1  meter  from  the 
sphere.  This  was  done  by  using  the  following  equation: 

V  xr 

Normalized  Scattering  =  - 

^incident 

where  r  is  the  distance  between  the  sphere  and  the  receiver  and  is  given  by 


V  =V'  X 

incident  incident 


R  +  0.15 
R 


where  is  the  measured  incident  signal  (approximately  15  cm  behind  where  the 

sphere  actually  was)  and  R  is  the  transmitter  to  sphere  distance. 


B.  THEORETICAL  DATA 

The  values  for  the  theoretical  data  were  obtained  from  two  FORTRAN  programs 
provided  by  Kargl.  One  program  calculates  the  theoretical  monostatic  data  results  and  the 
other  the  bistatic  results  for  a  saturated  poro-elastic  sphere  in  a  saturated  poro-elastic 
medium. 

The  program  requires  almost  30  inputs.  These  inputs  are  the  material  properties  of 
the  external  fluid  and  poro-elastic  medium,  both  water  in  this  experiment,  and  the  internal 
fluid  (water)  and  poro-elastic  medium  of  the  scatterer.  The  internal  poro-elastic  medium 
was  either  porous  glass  or  the  aluminum.  Inputs  to  the  programs  can  be  found  in 
Appendix  C. 
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c.  COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  DATA 


1.  Scattering  from  the  Aluminum  Sphere 

Figure  8  shows  a  comparison  of  the  theoretical  and  experimental  data  for 
monostatic  scattering  from  the  aluminum  sphere.  Experimental  data  points  are  marked  by 
X  and  connected  with  a  dashed  line.  The  dashes  are  only  to  guide  the  eye  and  are  not  an 
attempt  at  a  curve  fit.  After  the  experimental  data  were  taken  the  sphere  was  weighed  and 
measured  to  verify  the  properties  for  input  into  Kargl’s  program.  It  was  found  that  the 
aluminum  sphere  weighed  approximately  30%  more  than  it  should.  The  sphere  was  not  a 
solid  aluminum  sphere,  it  was  an  aluminum  shell  covering  a  core  of  some  heavier 
unknown  material.  The  main  features  in  the  experimental  data  agree  very  well  with  the 
theoretical  data  below  90  kHz  except  that  the  experimental  values  are  shifted  to  the  right 
by  5-10  kHz.  Beyond  90  kHz  the  features  match  in  relative  position  but  the  experimental 
data  has  a  lower  amplitude.  Figure  9  shows  the  data  adjusted  so  that  it  spans  20  to  150 
kHz.  This  shows  a  much  better  match  between  the  theoretical  and  experimental  values  at 
frequencies  below  90  kHz.  Since  the  sphere  was  determined  not  to  be  solid  aluminum 
these  results  are  considered  to  be  in  reasonable  agreement  with  the  predicted  values  and 
show  that  the  experimental  procedure  is  valid. 

For  the  bistatic  case,  shown  in  Figure  10  and  Figure  1 1,  the  structure  and 
amplitudes  do  not  agree  well.  Again,  this  is  not  unexpected  since  the  sphere  was  not  solid 
aluminum.  However,  because  the  lobing  is  clearly  defined  in  the  experimental  data,  it 
was  again  considered  that  the  experimental  procedure  was  valid. 
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Figure  8.  Normalized  monostatic  scattering  amplitudes  from  the  aluminum  based  sphere. 


Figure  9.  Normalized  monostatic  scattering  amplitudes  from  the  aluminum  based  sphere  with 

data  ‘stretched’  between  20  and  150  kHz. 


Figure  10.  Normalized  bistatic  scatt( 


Figure  11,  Normalized  bistatic 


Figure  12.  Normalized  monostatic  scattering  amplitudes  from  the  100  pm  sphere. 

2.  Scattering  from  the  100  p,m  Porous  Glass  Sphere 

Figure  12  shows  the  results  for  monostatic  scattering  from  the  100  jim  sphere. 
Most  of  the  features  in  the  experimental  data  can  be  found  to  correspond  fairly  well  to 
features  in  the  theoretical  model  up  to  about  80  kHz  and  less  so  above  80  kHz.  Note  that 
the  amplitude  of  the  experimental  data  is  lower  than  the  theoretical  data.  This  agrees  with 
the  bistatic  data,  which  consistently  shows  that  the  amplitudes  of  the  main  lobes  are 
lower  than  the  theoretical  prediction  even  when  the  side  lobes  have  similar  amplitudes  as 
shown  in  Figures  13  through  17. 

Figures  13  through  17  show  the  theoretical  and  experimental  bistatic  scattering 
data  sets.  There  is  good  agreement  for  the  30,  60  and  90  kHz  data  sets  except  for  the 
amplitudes  of  the  main  lobes.  The  120  kHz  data  matches  well  and  the  150  kHz  data  only 
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agrees  in  overall  magnitude  of  the  scattering.  This  is  believed  to  be  due  to  the  increased 
sensitivity  of  the  beam  pattern  to  the  material  properties.  In  varying  the  values  used  as 
input  for  Kargl’s  program,  it  was  noted  that  the  beam  pattern  varied  slightly  at  30  kHz 
and  increased  in  variation  as  frequency  was  increased.  The  resulting  beam  pattern  varied 
significantly  at  150  kHz.  This  is  the  most  likely  cause  of  the  disagreement. 


0° 


Figure  13.  Normalized  bistatic  scattering  at  30  kHz  from  100  gm  sphere. 
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Figure  16.  Normalized  bistatic  scattering  at  120  kHz  from  100  ^m  sphere. 
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Figure  17.  Normalized  bistatic  scattering  at  150  kHz  from  100  pm  sphere. 
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Figure  18.  Normalized  monostatic  scattering  from  the  500  |im  sphere. 

3.  Scattering  from  the  500  |im  Porous  Glass  Sphere 

Monostatic  scattering  for  the  500  |im  sphere  is  shown  in  Figure  18.  It  is  possible 
to  see  some  correspondence  between  features  below  90  kHz  but  above  this  frequency  no 
match  can  be  found. 

Figures  19  through  23  show  the  bistatic  results  for  the  500  |im  sphere.  There  is  no 
correspondence  between  the  experimental  and  theoretical  values.  It  is  believed  that  this  is 
due  to  defects  in  the  500  pm  sphere.  Prior  to  the  experiment  it  was  noted  that  this  sphere 
has  a  0.9  x  0.5  x  0.25  centimeter  gouge  in  its  surface.  The  surface  is  also  less  consistent 
in  texture  from  one  area  to  another.  Also  the  surface  is  more  susceptible  to  crumbling 
than  the  100  pm  sphere.  Finally,  both  spheres  were  weighed  and  measured  and  their 
densities  calculated.  Both  were  found  to  have  a  density  of  about  1550  kg/m\  This  density 
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corresponds  with  the  100  jim  bar  used  by  LT,  Huskey  for  determining  the  material 
properties  of  the  100  |im  sphere.  However,  the  density  of  the  500  pm  bar  used  by  LT. 
Huskey  was  1513  kg/m^  (Huskey,  1993)  which  is  about  2.5%  lower  than  for  the  500  pm 
sphere. 

Based  on  the  material  condition  of  the  500  pm  sphere  and  on  the  good  results 
obtained  with  the  100  pm  sphere  it  is  believed  that  the  500  pm  sphere  is  defective  and 
that  this  caused  the  disagreement  between  the  theoretical  and  experimental  data. 


0° 


Figure  19.  Normalized  bistatic  scattering  at  30  kHz  from  500  pm  sphere. 
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Figure  20.  Normalized  bistatic  scattering  at  60  kHz  from  500  p,m  sphere. 
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Figure  21.  Normalized  bistatic  scattering  at  90  kHz  from  500  u,m  sphere. 
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V.  CONCLUSION  AND  RECOMMENDATIONS 


A.  CONCLUSION 

The  scattering  of  sound  from  an  aluminum  based  sphere  and  from  porous  glass 
spheres  composed  of  100  |xm  and  500  |im  diameter  glass  beads  was  measured. 

Monostatic  and  bistatic  measurements  were  made  for  backscatter  at  30,  60,  90,  120  and 
150  kHz.  Experimental  data  were  compared  with  theoretical  data  computed  by  a 
FORTRAN  program  written  by  Kargl.  The  program  was  based  on  the  theory  developed 
by  Kargl  and  Lim  for  the  scattering  of  sound  from  a  fluid  saturated  poro-elastic  sphere  in 
a  saturated  poro-elastic  medium. 

Material  properties  of  the  spheres  were  measured  for  input  into  Kargl’ s  program. 
The  elastic  moduli  of  a  cylindrical  bar  composed  of  300  |im  diameter  bonded  glass  beads 
were  measured  and  used  to  determine  the  bar’s  Poisson’s  Ratio.  This  Poison’s  ratio  was 
used  to  determine  the  value  of  the  bulk  moduli  of  the  100  |xm  and  500  pm  porous  glass 
spheres.  Other  material  properties  of  these  spheres  had  been  previously  measured  by  LT. 
Huskey. 

Comparison  of  the  experimental  data  to  the  theoretical  from  Kargl’ s  program 
yielded  reasonable  results  for  the  aluminum  based  sphere.  Very  good  results  were 
obtained  for  the  100  pm  sphere.  The  measurements  at  the  lower  frequencies  agreed  more 
closely  than  those  at  the  higher  frequencies.  This  is  due  to  the  increased  sensitivity  of  the 
spatial  structure  of  the  scattering  at  higher  frequencies  to  slight  variations  in  the  material 
properties  of  the  spheres.  The  main  lobe  was  found  to  be  lower  in  amplitude  than 
predicted  at  all  frequencies.  No  explanation  of  this  could  be  found  at  this  time.  It  is 
speculated  that  this  may  be  the  result  of  slight  inhomogeneity  in  the  composition  of  the 
sphere.  At  higher  frequencies  it  may  be  due  to  error  in  the  axial  alignment  of  the  source, 
sphere,  and  receiver.  The  500  pm  sphere  produced  extremely  poor  results  and  it  is 
believed  that  this  sphere  is  defective,  specifically  it  is  thought  to  be  non-homogene ous. 
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B.  RECOMMENDATIONS 


Scattering  measurements  at  lower  frequencies  is  recommended.  This  would 
require  a  larger  tank  than  is  present  at  the  Naval  Postgraduate  School.  This  would  allow 
measurements  to  be  made  where  the  fluid’s  viscosity  would  affect  the  flow  of  the  fluid 
through  the  pores  in  the  spheres. 

For  measurements  at  higher  frequencies  the  material  properties  of  the  spheres 
need  to  be  known  to  fairly  high  accuracy.  For  this  to  be  accomplished  new  spheres  should 
be  made  at  the  same  as  the  time  cylindrical  samples  are  made.  The  cylindrical  samples 
should  have  a  length  to  diameter  ratio  greater  than  15:1  to  allow  accurate  measurement  of 
the  elastic  moduli.  The  surfaces  of  the  cylinders  should  also  be  machined  so  that  they  are 
consistent  with  the  sphere’s  surfaces.  Also,  for  measurements  at  higher  frequencies,  a 
more  accurate  method  for  positioning  the  receiver  is  needed.  Measurements  should  be 
taken  at  sufficiently  small  intervals  to  adequately  define  the  structure  of  the  side  lobes. 

Finally,  an  investigation  into  the  non-uniform  variation  of  the  scatterer’s  material 
properties  should  be  performed.  This  would  determine  their  affect  on  the  structure  of  the 
beam  pattern  from  the  sphere  which  may  explain  the  poor  results  obtained  with  the 
500  p.m  sphere. 
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APPENDIX  A.  MONOSTATIC  DATA  RESULTS 
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Frequency 
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3,030.4 

0.0054 

119.62 

6,078.0 

0.0141 

BcairJW 

0.0137 

18.33 

3,179.8 

0.0041 

142.32 

6,240.0 

0.0163 

mnsm 

0.0161 

32.47 

3,251.4 

0.0071 

141.69 

6,258.3 

0.0162 

130 

78.00 

3215.0 

28.93 

3,144.4 

0.0066 

122.84 

6,090.8 

0.0144 

132 

70.86 

3032.4 

0.0168 

— 

2,935.0 

— 

119.41  ■ 

5,806.3 

KBESH 

134 

58.31 

2902.6 

0.0145 

16.47 

2,787.2 

0.0042 

126.25 

5,488.0 

IKI^I 

136 

44.83 

2781.8 

0.0116 

20.30 

2,653.4 

0.0055 

120.71 

5,142.0 

138 

38.13 

■IE 

0.0046 

102.76 

4,747.1 

140 

■eb 

■EB 

0.0045 

81.74 

4,352.2 

142 

0.0135 

18.04 

0.0065 

66.25 

3,965.6 

144 

33.93 

1943.0 

0.0126 

18.07 

0.0072 

55.41 

3,611.2 

0.0110 

146 

1796.5 

0.0108 

148 

1655.3 

0.0087 

mmm 

■Bl 

KB 

3,063.8 

■EB 

150 

1481.7 

0.0065 

6.84 

1,417.4 

0.0035 

LiHi- 

2,757.6 

0.0109  1 
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appendix  b.  bistatic  daia  results 


Aluminum 


2 


Incident 


Scattered 

Amplitude 

(mV) 

Normalized 

Amplitude 

1.84 

2.76 

0.0266 

2.96 

0.0286 

2.99 

0.0288 

2.98 

0.0288 

US 

Bl^l 

hesbi 

1.65 

0.0159 

1,72 

0.0166 

1.57 

0.0152 

1.79 

0.0173 

2.15 

0.0208 

2,05 

1.63 

Scattered 

Amplitude 

(mV) 


0.67 

0.69 


Normalized 

Amplitude 


0.0066 

0.0067 

0.0138 


500  \im 

Scattered 

Amplitude 

(mV) 

Normalized 

Amplitude 

1.40 

1.11 

0.64 

0.86 

0.95 

1.50 

L  0.0146 

0.73 

0.0071 

1.07 

0.0104 

0.92 

0.0090 

1.13 

1.22  1 

0.34  1 

0.25 

0.0024 

1.41 

0.0137 

1.26  1 

0.0123 

1.65 

0.0161 

1.73 

0.0168 

2.38 

0.0232 

2.94 

0.0287 

4.01 

0.0391 

3.61 

0.0351 

3.64 

0.0354 

1.33 

0.0129 

1.74 

0.0169 

1,41 

0.0137 

1.46 

0.0142 

1.43 

0.0139 

1.60 

0.0156 

1.31 

0.0128 

1.46 

0.0142 

1.79 

0.0173 

1.92 

0.0185 

3.42 

0.0329 

2.59 

1,33 

74.22 

1.80 

0.0175 

1.56 

0.0151 

2.27 

0.0221 

2.08 

0.0202 

1.71 

0.0167 

1.39 

0.0136 

0.55 

0.0053 

0.15 

0.0015 

0.67 

0.0066 

1,45 

0.0142 

2.15 

0.0209 

2.34 

0.0228 

2.12 

0.0207 

.21 

0.0118 

.40 

0.0137 

.26 

0.0122 

.51 

0.0147 

.08 

0.0106 

.53 

0.0052 

0.17 

0.0017 

1.15 

0.0112 

1.18 

0.0114 

1.39 

0.0136 

1.35 

0.0131 

1.00 

0.0098 

0.78 

0.0076 

0.91 

0.0089 

73.51 
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60  kHz 


34 


90  kHz 


35 


120  kHz 


230 


235 


240 


245 


250 


255 


260 


265 


270 


Incident 


100  fim 

Scattered 

Amplitude 

(mV) 

Normalized 

Amplitude 

33.20 

39.00 

29.31 

33.11 

23.39 

Scattered 

Amplitude 

(mV) 


3877 


34.00 

0.0077 

35.89 

0.0082 

18.77 

0.0043 

17.43 

1  0.0040 

43.54 

0.0099 

54.23 

0.0123 

53.13 

0.0121 

17.51 

0.0040 

41.77 

0.0095 

44.56 

0.0101 

31.83 

0.0072 

32.51 

0.0074 

40.28 

0.0092 

36.76 

0.0084 

18.17 

0.0041 

2 

0.0105 

3 

0.0063 

Normalized 

Amplitude 


0.0106 


28.49 

11.18 

20.21 

0.0055 

29.50 

0.0081 

26.58 

0.0073 

36.97 

0.0101 

28.56 

0.0078 

17.68 

0.0048 

19.39 

0.0053 

23.07 

0.0063 

24.95 

0.0068 

18.48 

0.0051 

24.74 

0.0068 

27.82 

0.0076 

26.91 

0.0074 

i  0.0083 

0.0072 

0.0087 

22.21 

0.0061 

41.32 

45.13 

26.00 

0.0071 

36.59 

0.0100 

35.17  i 

0.0096 

29.53 

0.0081 

31.48 

0.0086 

26.86 

0.0074 

27.25 

0.0075 

0.0086 

IEOH 

0.0077 

1  39.94 

0.0109 

mtimm 

1  16.91 

0.0046 

0.0077 

0.0117 

1  2612.0 

150  kHz 


Aluminum 

100  |Lim 

Scattered 

Scattered 

Amplitude  Normalized 

Amplitude 

Normalized 

Angle 

(mV)  Amplitude 

(mV) 

Amplitude 

260 


265 


270 


Incident 


500  |im 

Scattered 

Amplitude 

(mV) 

Normalized 

Amplitude 

11.00 

0.0135 

15.08 

0.0185 

13.85 

0.0170 

12.54 

0.0154 

27.40 

0.0336 

0.0026 

0.0045 

7.66 

0.0037 

9.19 

0.0044 

15.04 


30.90 


4.89  I 


23.90 


0.0045 

0.0037 

13.59 

0.0066 

15.24 

0.0074 

0.23 

0.0001 

4.73 

0.0023 

11.87 

0.0057 

8.62 

0.0042 

0.0074 

17.47 

0.0214 

33.76 

0.0413 

18.08 

0.0221 

0.0242 

32.76 

8.70 

0.0107 

16.79 

0.0206 

8.93 

0.0109 

33.38 

0.0409 

1.36 

0.0017 

18.24 

0.0223 

9.89 

0.0121 

8.46 

0.0104 

CVJ 

cd 

0.0075 

21.72 

0.0266 

25.74 

15.17 

0.0186 

545.3J 

571.48 

16.71 

0.0081 

19.46 

0.0094 

11.70 

3.70 

5.96 

4.63 

6.58 

13.05 

6.60 

15.67 

6.94 

0.0034 

9.36 

0.0045 

11.79 

10.84 

6.36 

1 

15.67 

21.24 

16.22 

1.75 

■ 

8.85 

0.0047 

6.76 

0.0036 

5.49 

4.83 

9.90 

Wi 


9.17 

0.0049 

12.60  ! 

0.0067 

iIQ] 

Wi 


9.24 

13.97 

13.83 

6.39 

'  15.89 

;  1337.5 
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APPENDIX  C.  INPUTS  TO  KARGL’S  PROGRAM 


i  Monostatic  Inputs  | 

Aluminum 

100  |im 

500  jim 

External  Fluid  (Water) 

Density,  p, 

998.665 

998.665 

998.665 

Bulk  Modulus,  K, 

2.17293x10® 

2.17293x10® 

2.17293x10® 

Viscosity,  Tj 

0.0 

0.0 

0.0 

Internal  Fluid 

Density,  p„ 

2700 

998.665 

998.665 

Bulk  Modulus,  Kfo 

8.078x10^° 

2.17293x10® 

2.17293x10® 

Viscosity,  rjo 

0.001 

0.001 

0.001 

External  Medium  (Water) 

Density,  Pf 

998.665 

998.665 

998.665 

Solid  Bulk  Modulus, 

(2.17293x10®,  0.0) 

(2.17293x10®,  0.0) 

(2.17293x10®,  0.0) 

Lattice  Bulk  Modulus,  Kg 

(2.25x10®,  0.0) 

(2.17293x10®,  0.0) 

(2.17293x10®,  0.0) 

Shear  Modulus,  p 

(1.0,  0.0) 

(1.0,  0.0) 

(1.0,  0.0) 

Tortuosity,  a 

1.65 

1.65 

1.65 

Porosity,  p 

0.999999 

0.999999 

0.999999 

Permeability,  /c^ 

1.0 

1.0 

1.0 

Internal  Medium  (Sphere) 

Density,  po 

2700 

2231 

2231 

Solid  Bulk  Modulus,  K^o 

(8.078x10'°,  0.0) 

(3.5x1 0'°,  0.0) 

(3.5x1 0'°,  0.0) 

Lattice  Bulk  Modulus,  Kgo 

(8.078x10'°,  0.0) 

(2.96x10®,  0.0) 

(2.87x10®,  0.0) 

Shear  Modulus,  Po 

(2.677x10'°,  0.0) 

(2.81x10®,  0.0) 

(2.72x10®,  0.0) 

Tortuosity,  ag 

1.65 

1.65 

1.65 

Porosity,  po 

0.321 

0.306 

0.305 

Permeability, 

1.0x10-'° 

6.53x10-'® 

5.74x10-" 

Miscellaneous 

ap 

1.0x10-® 

1.0x10-® 

1.0x10-® 

Spo 

1.0x10® 

1.0x10-® 

1.0x10-® 

dfreq 

100 

100 

100 

max.  freq. 

150000 

150000 

150000 

radius 

0.0382 

0.0344 

0.0346 

distance 

0.75 

0.75 

0.75 

hash 

0 

0 

0 

nstart 

0 

0 

0 

nend 

75 

75 

75 
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Bistatic  Inputs 


External  Fluid  (Water) 

Density,  pf _ 

Bulk  Modulus.  K, _ 

Viscosity,  rj  _ 

Internal  Fluid  _ 

Density,  pm  _ 

Bulk  Modulus,  _ 

Viscosity,  rio _ 

External  Medium  (Water) 
Density,  p, 

Solid  Bulk  Modulus, 

Shear  Modulus,  _ 

Porosity,  p  _ 

Permeability,  _ 

Internal  Medium  (Sphere) 

Density,  pp _ 

Solid  Bulk  Modulus,  Ksp 
Shear  Modulus,  tip 

Porosity,  Po _ _ 

Permeability,  _ 


Miscellaneous  _ 


nfreq  _ _ 

freqmin _ _ 

dfreq 


Aluminum 

998.6^ 

2.17293x10^ 

0^ 

2700 

8.078x10'° 

1.0 _ 

998.665 

(2.17293x10°,  0.0) 
(1.0,  0.0) 
0.999999 


2700 

(8.078x10'°,  0.0) 
(2.677x10'°,  0.0) 
1.0x10° 
1.0x10-'° 


30000 

120000 


100  |um 


998.665 

2.17293x10° 


998.665 

2.17293x10° 

0.001 

998.665 

(2.17293x10°,  0.0) 
(1.0,  0.0) 
0.999999 


(3.5x1 0'°,  0.0) 
(2.6x1 0'°,  0.0) 
0.306 
6.53x10-'° 


30000 

30000 


500  |im _ 

998.665 

2.17293x10° 

0.0 

998.665 

2.17293x10° 

0.001 

998.665 

(2.17293x10°,  0.0) 
(1.0,  0.0) 
0.999999 


(3.5x1 0'°,  0.0) 
(2.6x1 0'°,  0.0) 
0.305 
5.74x10-" 


30000 

30000 


distance 

hash 


nstart 


#  of  angles _ 

Experimental  K^p 
Experimental  pp 


(2.96x10°,  0.0) 
(2.81x10°,  0.0) 


(2.87x10°,  0.0) 
(2.72x10°,  0.0) 
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